Noname manuscript No. 

(will be inserted by the editor) 



Hermann Freire • Eberth Correa 

A functional generalization of the 
field-theoretical renormalization group 
approach for the single-impurity 
Anderson model 



Received: date / Accepted: date 



Abstract We apply a functional implementation of the field-theoretical renor- 
malization group (RG) method up to two loops to the single-impurity Ander- 
son model. To achieve this, we follow a RG strategy similar to that proposed 
by Vojta et al. [Phys. Rev. Lett. 85, 4940 (2000)], which consists of defin- 
ing a soft ultraviolet regulator in the space of Matsubara frequencies for 
the renormalized Green's function. Then we proceed to derive analytically 
and solve numerically integro-differential flow equations for the effective cou- 
plings and the quasiparticle weight of the present model, which fully treat 
the interplay of particle-particle and particle-hole parquet diagrams and the 
effect of the two-loop self-energy feedback into them. We show that our re- 
sults correctly reproduce accurate numerical renormalization group data for 
weak to slightly moderate interactions. These results are in excellent agree- 
ment with other functional Wilsonian RG works available in the literature. 
Since the field-theoretical RG method turns out to be easier to implement at 
higher loops than the Wilsonian approach, higher-order calculations within 
the present approach could improve further the results for this model at 
stronger couplings. We argue that the present RG scheme could thus offer 
a possible alternative to other functional RG methods to describe electronic 
correlations within this model. 
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1 Introduction 

The single-impurity Anderson model (SIAM) has always been a central model 
in the field of strongly correlated systems pQ. It was originally proposed to 
describe microscopically the effect of a localized impurity embedded in a 
metallic host on its conduction electrons. Recently, the importance of this 
model has increased even further. Thanks to ground-breaking advances in 
nanotechnology, it is now possible to simulate experimentally this model in 
nanostructure devices such as quantum dots coupled to metallic leads [2J. 
This has opened a new arena for exploration of correlation effects in many- 
body systems. On the theoretical side, this renewed interest in the model is 
also related to the discovery that the paradigmatic Hubbard model in the 
appropriate limit of infinite spatial dimensions [3] (or infinite coordination 
number) can in fact be mapped onto an effective SIAM with an additional 
self-consistency condition pQ. This discovery led to the development of the 
Dynamical Mean-Field Theory (DMFT) which, combined with first-principle 
methods such as the Density Functional Theory [5] to determine the elec- 
tronic structure, is currently being applied to realistic strongly correlated 
materials with good success. 

However, since it is difficult to solve the SIAM in its most general form 
(except for some special conditions where it can be solved with the help 
of the Bethe ansatz [5] or the Numerical Renormalization Group (NRG) 
[7]), this represents nowadays an obstacle to some extent for the systematic 
application of DMFT to many strongly-correlated materials. There is a lack 
of reliable methods to study some of the complicated generalizations of the 
SIAM generated by DMFT calculations, especially in the strong-coupling 
regime. To tackle this problem, it is important to set up new theoretical 
schemes which offer the advantage of being flexible in the implementation or 
generalization and which require preferably low computational effort. In this 
respect, we mention that a major step forward has been recently achieved 
via diagrammatic quantum Monte Carlo which turns out to be free from the 
fermionic sign problem [5]. 

Another promising semi-analytical method developed in recent years is 
based on the functional generalization of the renormalization group (RG) 
approach. Several different implementations of this technique are already 
available in the literature, both in the Wilsonian version [9, f01[TTlll2lll4l[T51 
UMIMrn and the field-theoretical version [TSlEIlinil^J . Th ey have been 
applied to many electronic models ranging from the two-dimensional (2D) 
Hubbard model to quantum impurity models such as the SIAM or the closely 
related Kondo model. Concerning the 2D Hubbard model, for instance, some 
of these functional RG works successfully reproduced within a static ap- 
proximation both at one-loop [THll5llT51rT9] and two-loop |23[|2~I] levels an 
antiferromagnetic phase in the model near half-filling and the onset of a 
d 2; 2_j ( 2-wave singlet superconducting phase away from half-filling for weak- 
to-moderate couplings. This agrees qualitatively with the physics displayed 
by the cuprate superconductors and gives further support to the point of 
view that this 2D model might indeed capture some important features of 
these strongly-correlated materials. 
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Given these encouraging results, several functional RG studies have also 
been initiated nowadays to discuss quantum impurity problems [2"5"ll2"o] such 
as the SIAM with the important goal of reproducing some key aspects of this 
model, most notably, its dynamical properties for all interaction strengths 
and also the correct emergence of the exponential Kondo scale at strong 
coupling. In addition to providing useful benchmarks to the functional RG 
scheme - more specifically, to the most common approximations used in these 
approaches - these studies could also provide alternative and computationally 
amenable semi-analytical quantum impurity solvers which, as we have seen 
before, are a central ingredient in the context of DMFT calculations. 

In this respect, Karrasch et al. implemented a Wilsonian functional RG 
study at one-loop level [25] of several important quantities of the model such 
as the effective mass, the static spin susceptibility and the spectral func- 
tion. They concluded that while their approach yielded very good results for 
weak couplings, they still could not reproduce the exponentially small Kondo 
asymptotics at higher couplings. Besides that, due to difficulties (ill-posed 
problem) with the stability of the analytic continuation of their functional 
RG data for the the single-particle Green's function, they were not able to 
establish conclusively the formation of the so-called Hubbard satellite bands 
which emerge at high energies in the impurity spectral function for strong 
coupling. Recently, two different works reported some progress on these prob- 
lems: Using a Hubbard- Stratonovich decoupling of the exact Fermi-Bose RG 
equations, Isidori et al. managed to obtain a three-peak structure in the spec- 
tral function of the model at finite temperatures [27] in qualitative agreement 
with precise NRG calculations for intermediate couplings. Also, Jakobs et al. 
applying a functional RG scheme within the Keldysh formalism }28] obtained 
qualitative similar results for the SIAM in the equilibrium situation for weak- 
to-moderate interaction. Despite that, both approaches could not obtain an 
exponential narrowing of the quasiparticle peak at strong coupling, which 
suggests that the Kondo scale emergence is still not completely included at 
the level of approximation used in those works. In this sense, a next logical 
step would be to check if adding higher-order terms would improve further 
these results for the model. It is interesting to mention that higher-order 
terms are also recognized as an important contribution in order to compute 
analytic results for transport quantities and dynamical correlation functions 
of the Kondo model out of equilibrium within the so-called real-time renor- 
malization group method in frequency space I291I301I5T] . 

In the present work, we implement a two-loop calculation of the the self- 
energy together with a calculation of the the renormalized coupling functions 
for the SIAM at zero temperature using a functional generalization of the 
standard field-theoretical RG approach. This analysis is important in view of 
the fact that it includes self-energy effects in all coupled integro-differential 
RG flow equations. In addition to this, it also includes, by construction, 
the dynamics of all frequency-dependent quantities in the model. Since the 
functional field-theoretical RG is not widely known among condensed matter 
physicists we present the method at length and in full detail. As an initial 
test for this RG scheme, we present both analytical and numerical results for 
the renormalization of the effective coupling functions and the quasiparticle 
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weight in the low-energy limit for all interaction strengths. Afterwards, we 
proceed to benchmark our results against exact (or highly accurate) data for 
this model such as the Bethe ansatz and NRG. 

This paper is structured as follows. In Sec. II, the single impurity An- 
derson model that we want to investigate is introduced. Next (in Sec. Ill), 
we explain the functional field-theoretical RG methodology used to study 
the model. The two-loop RG flow equations are derived explicitly in Sec. 
IV. We then proceed to solve numerically the coupled integro-differential RG 
equations and discuss the main results in Sec. V. Finally, we present our 
conclusions from this study. 



2 Model 

The Hamiltonian of the SIAM is given by 



where Ck CT annihilates a noninteracting conduction electron with momentum 
k, energy dispersion et, and spin projection a, whereas d a annihilates a 
localized electron with energy Ed and on-site repulsion U. The hybridization 
between the localized electrons and the conduction electrons is characterized 
by the amplitude Vk- 

If we use a coherent-state functional integral representation of Eq. ([1]) at 
constant chemical potential fi and integrate out the conduction electrons, the 
corresponding model at T — becomes readily described by the action 



where £o = Ed— \i and J = J In addition, d a and d a are, respectively, the 
creation and annihilation Grassmann fields for localized electrons with spin 
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projection a and with the hybridization function given by A(iu)) — ■ 
The above action then defines our bare quantum field theory which must be 
appropriately regularized in the UV regime (more details on the regulariza- 
tion procedure in the next section). For simplicity, we shall consider here 
the SIAM in the so-called wide-band limit such that A(iui) — — iAsga(oj), 
where A is a constant. It is also convenient to make the following choice of 
variables, i.e. £o = —U/2. When this is done, the model becomes completely 
particle-hole symmetric. 




+ X^k4c kff + v k 4«A) + ud\d t d\d h 



(1) 




(2) 
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If we treat the interaction of Eq. ([2]) within the self-consistent Hartree- 
Fock approximation, we obtain the renormalized excitation energy = 
£o + S hf (uj), where 

S BF (u) = U f G hf {ilu-ilo'), (3) 

with GHF(i<j->) = l/\iuj — + iAsgn(u)} being the self-consistent Hartrec- 
Fock impurity Green's function of the model. Since, for frequency-independent 
interaction, one obtains that Ehf(uj) = U/2, this latter contribution pre- 
cisely cancels £o hi the particle-hole symmetric case. However, as will become 
clear soon, despite the fact that the initial interaction of the SIAM is indeed 
a constant, it will acquire a frequency-dependent structure under the RG 
transformation. For this reason, in order to preserve the particle-hole sym- 
metry of the model during the whole RG flow, we must go back to Eq. @ 
and define a frequency-dependent impurity energy £o = ^o(^), such that the 
condition £o( w ) + /„/ U(uj,uj',uj')GHF(iw — iw') = is always enforced for 
this case. As a consequence of this, the corresponding Hartree-Fock impurity 
propagator will naturally take the standard result 

Ghf(lo) = ^— --. (4) 

iui + iZAsgn(w) 



3 RG Methodology 

The outline of the field-theoretical RG method can be generally summa- 
rized in the following way (for more details, see, e.g., Ref. [32]). Typically, in 
correlated systems if one applies a naive perturbation theory for a given elec- 
tronic model, divergences or non-analyticities often appear in the low-energy 
limit at the calculation of several physical quantities of the model such as, 
e.g., irreducible two-particle vertices, the self-energy and also several order- 
parameter susceptibilities. This result usually implies that the perturbation 
theory setup is not appropriately formulated. The field-theoretical RG ap- 
proach circumvents this problem by means of a reorganization of the terms 
in the conventional perturbative approach. More precisely, by rewriting the 
experimentally unobserved bare quantities of the microscopic field theory 
model (such as, e.g., the fermionic fields and the coupling parameters) in 
terms of the corresponding renormalized parameters, it is possible to con- 
struct a new renormalized perturbation theory which is, by construction, 
convergent in the low-energy limit. The bare quantities are always defined at 
the ultraviolet (UV) scale of the model, whereas the renormalized parameters 
are conveniently defined at a physically measurable low-energy scale. Since 
the renormalized parameters are the ones that are physically probed from 
the experimental point of view, they must come out necessarily UV cutoff- 
independent. This regularization procedure has to be implemented order by 
order in perturbation theory. If this program is successfully accomplished, 
then the field theory model is said to be properly renormalized. 

It is interesting to note however that, unlike the situation for many lat- 
tice models such as, e.g, the 2D Hubbard model, where divergences and 
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Fig. 1 The self-energy diagram of the SIAM at two-loop order. The thick lines 
represent the interacting Green's function of the model. 



non-analyticities in Feynman diagrams often abound, perturbation theory 
is regular for the SIAM. Therefore, one might naively think, at first sight, 
that there should be no intrinsic need in this model for a RG-based theo- 
retical attack. This view is overly restrictive though, especially due to the 
fact that the emergence of an exponentially dependent (in U) Kondo scale 
at the model indeed motivates many-body resummations of certain classes 
of diagrams. The RG framework is an ideal tool for achieving this goal in an 
unbiased manner. 

In this work, we implement a functional field-theoretical RG approach up 
to two-loops for the SIAM in order to calculate the dynamics of all frequency- 
dependent renormalized quantities of the model. This analysis is important 
for the following reasons: Firstly, despite the fact that the frequency depen- 
dencies of the renormalized parameters could be classified as irrelevant in a 
power counting sense, they do have a crucial effect on the low-energy be- 
havior of the model especially at stronger couplings. Secondly, in order to 
evaluate experimentally relevant dynamical quantities such as the spectral 
density function and the dynamical susceptibilities of the model, it is essen- 
tial to include such a generalization from the outset in all physical parameters 
of the RG approach. Lastly, we briefly mention that finite-energy properties 
of the model are also important if one is willing to extend the functional 
field-theoretical RG approach to the non-equilibrium situation, e.g., to cal- 
culate transport properties in quantum dots (see, for instance, Ref. |28j for 
a Wilsonian RG approach to this problem). 

To begin our analysis, we describe the underlying strategy of our RG 
scheme. We shall follow a RG method that has some similarities with the 
approach proposed by Vojta et al. in a different context [33] . namely to 
discuss possible critical points in high-T c superconductors. In our case, in 
order to implement a functional generalization of the field-theoretical RG 
method and also to avoid possible artificial non-analyticities generated by 
the regularization procedure in the SIAM, it is better not to impose a sharp 
UV cutoff in the field-theory model defined by Eq. ((2J. For this reason, we 
define a soft regulator if (|u;|/yl) in the space of Matsubara frequencies for 
the fully interacting Green's function (where A stands for the UV cutoff of 
the model) . The exact impurity propagator should then be given by 

G a (ilo) = - 1 K{\u\/A), (5) 

ioj + iAsgn{uj) — Zyyi(oj) 

where Sa{^>) = Ea{u>) — ^a_ F i^) 1S * ne self-energy of the SIAM without the 
Hartree-Fock term, since this latter contribution is already included in the 
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"noninteracting" propagator of the model as we have shown previously. In 
addition to this, K(y) must be some smooth decaying function which must 
satisfy K(0) = 1. Following Ref. [53], we conveniently choose K(y) = e~ v . 
This choice of the UV regulator turns out also to be more amenable from 
both analytical and numerical point of view. (We parenthetically note here 
that a sharp cutoff would correspond to a Heaviside step function K(y) = 
9(y — 1).) As we will discuss shortly, to derive the RG flow equations, in 
which degrees of freedom with Matsubara frequencies lying between A and 
A — dA are successively integrated out, we must take a A (d/dA) derivative 
of all irreducible two-particle vertices and the self-energy of the model. 

Since we will be interested in the low-energy behavior of the SIAM - 
particularly in the possible emergence of the dynamically generated Kondo 
scale - we shall focus here on the low-frequency range \lu\ < A of the model. 
In this regime, it is reasonable to assume that the renormalized Green's 
function is given by the low-frequency Fermi liquid form 

Ga{^) = . m . 7 Za a . , K (MM) , (6) 

i.e. the self-energy of the model is given by S^iiuS) — iw{l — 1/Za), where 
Za = (1 — 917/1 (iw)/c?(iu;)| w= o) _1 is the quasiparticle weight. At weak- 
coupling, we know from second-order perturbation theory |34j that this quan- 
tity should be given by 

Z A ^o = l + ( 3 -x) i u /* A )* + -> ( 7 ) 

whereas in the strong-coupling limit U — > oo, from Bethe ansatz results [BJ, 
the quasiparticle weight is given asymptotically by 



Z A ^o ex P[~^ U / 8A ^ ■ ( 8 ) 

In what follows, we calculate explicitly the RG flow equations up to two- 
loop order for the SIAM. In order to do this, we choose a parametrization of 
the interaction parameter of the model in a way which manifestly incorpo- 
rates the energy and the spin conservation 



The RG methodology we apply now is standard and follows closely the 
field-theoretical method (see, e.g., Refs. [521155] ). We begin by first calculating 
the RG flow equation for the self-energy of the model at two-loop order as 
depicted in Fig. 1. As a result, we obtain 
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d A E A (iu) 



Ait 2 A 2 



+ 00 



+ OO 



du" 



Uf(w, OJ+Oj"- Lj',Lj")ljf(Lj', Lj" , LO + w" - J) 

+ 2U^(u),u) + u)" -lj',lj")U^(lj',lj",lj + lj" - lj') - Uf(u, uj" , u')Uf (oj + oj" - oj' , oj',oj) 

- 2Kft(w + oj" - u',u\ w)K^(w, oj", oj')] # (l^'IM) # (KIM) # (|w + w" - 

(|w'| + |w"| + |w + w" - u/|) 



[w' + Z A Asgn(oi')] [lj" + Z A Asgn(u")] [lj + lj" - lj' + Z A Asgn(Lj + lj" - lj')} 



(10) 



An important point we wish to emphasize here is that it is essential to use 
the interacting Green's function (i.e., which includes self-energy feedback) 
in the RG perturbative scheme instead of the Hartree-Fock Green's function 
in order to remove the unphysical Stoner instability in the model (for a 
discussion on this issue, see, e.g., Ref. [26]). If we now use the many-body 
definition displayed above for the quasiparticle weight factor Z A , we can 
straightforwardly derive the RG flow equation at two-loop order for this 
quantity as follows 

A d\nZ A 

A^T = v, (ii) 

where r\ = AZ A -^^ [d A E A {iojj\ u=Q . In this way, 77 is given by 



doj' I dLj" {[Uf(u,U + Lj" - Oj',Oj")Uf(uj' ,Lj" ,LJ + Lj" - Oj') 



Att 2 A 

+ 2U^(oj,oj+oj" -oj',oj")U^(oj',oj",oj + oj" -oj') - U({oj,oj",oj')U((oj+< 
j A . 



2Ul\(oj + lj" - lj 1 , lj', u)UA{w, lj", lj')] F a (lj, lj' , lj") ) 

) LJ—0 



- to 1 , u; , to) 
(12) 



where 



T a (lj,lj',lj") 



K{\lj'\/A)K{\lj"\/A)K{\lj + lj" - lj'\/A) 



(lj' + Z A Asgn(oj')){oj" + Z A A sgn(w"))(w + lj" - oj' + Z A Asgn(oj + lj" - lj')) 
d 

du 



„ ,. (\lj'\ + \lj"\ + \lj + lj" -lj'\) d ,. .. 

u -- lj" -oj'\- - — 1 1 ! ■ - — \lj + lj" - lj'\ 



A dLj 
\lj'\ + \lj"\ + \lj + lj" - w'|)(l + Z A Af- sgn(a; + lj" - lj')) 



(lj + lj" - lj' + Z A A sgn(w + lj" - lj')) 



(13) 



Next, we proceed to calculate the RG flow equations for the coupling 
functions of the SIAM. Since all the renormalized parameters are the phys- 
ically measurable quantities, they do not depend on the UV cutoff scale A. 
Therefore, we must set AdrR i i(oJx,oj2,oJz)/dA = for i =f and tl, where 
r^i(LJi,Lj2, 0J3) denote the renormalized irreducible two-particle vertices [llj . 
The corresponding Feynman diagrams for the coupling functions with two- 
loop self-energy feedback are schematically shown in Fig. 2. As a result, we 
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Fig. 2 The parquet diagrams for the renormalized irreducible two-particle ver- 
tices of the SIAM with two-loop self-energy feedback including the competition of 
both particle-particle and particle-hole contributions. The thick lines represent the 
interacting Green's function of the model. 

obtain the following coupled integro-differential RG equations which fully 
describe the interplay of both particle-particle and particle-hole vertex cor- 
rections in the model (the so-called parquet diagrams), i.e. 



uu* I I , ^2 J } AAA 

A — — — =A$(u}i,w 2 ,u} 3 )+Bf (wi,W2,w 3 ) + C£ (wi,^,^)) (14) 

and 



A — = J 4 ti (wi,w 2 ,w 3 ) + B n (u)i,u>2,u) 3 ) + C n (cji,w2,w3), 

(15) 

where LU4 =^1+^2—^3 and the functions A^-, B^, Aft , fift , and are 
given explicitly in Appendix A. The Eqs. (|TT|) . ([T4"|) . and (|15l) must be solved 
simultaneously in order to analyze the low-energy behavior of this model. 
Since it is impossible to obtain an analytical solution for these equations in 
this case, we must then resort to numerical methods. This will be done in 
the next section. 



4 Numerical Results 

We now devote this section to the numerical analysis of the integro-differential 
RG equations up to two-loop order derived above. Following previous works 
25,27], since we are mostly interested in the low-energy physics of the SIAM 
in order to reproduce the Kondo scale we discretize the Matsubara frequen- 
cies by introducing the geometric mesh 
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Fig. 3 (Color online) RG flow with two-loop self-energy corrections for some 
choices of external frequencies of the couplings U^(wi, ui-z, W3) and f-f (wi, W2, W3) 
for the SIAM (in units of A) as a function of the ratio A/ A. The bare coupling at 
the initial cutoff (Aq/A) = 10 was chosen to be (U/A) = 4. 



uj n = ±uj — n = 1...N, (16) 

a — 1 

where the free adjustable parameters ujq > and a > 1 determine the spac- 
ing of the discretized frequencies and N is the total number of frequencies to 
be used. For this discretization procedure, the low-frequency regime will be 
better resolved than the corresponding high-frequency sector, which is justi- 
fied since the Kondo peak emerges at lj = in the model. A typical choice 
of these parameters is, e.g., 

oj = 1(T 5 A a = 3.5, A = 10A, N = 30 points, 

which allows one to solve the flow equations at two-loop order with relatively 
low computational effort. Despite this moderate choice of number of points 
in the discretization procedure, we point out that our results show good 
convergence properties in the low-energy limit. This choice implies that we 
have to keep track of the RG flow of thousands of couplings altogether. Since 
the coupling functions in the RG equations need to be evaluated for arbitrary 
arguments - which in general do not coincide with a point of the mesh - an 
interpolation procedure has to be implemented. In this work, we shall use a 
linear interpolation. As a result of this, all the integrals in the flow equations 
will be solved by using standard Gaussian quadrature routines. Afterwards, 
the resulting set of coupled differential RG equations are then solved by 
means of high-precision fourth-order Runge-Kutta method. 

First, we focus our attention on the numerical solution of Eqs. (fT4)) and 
(fTS"]) for the renormalized couplings as a function of A. The choice of the initial 
conditions at A = Ao in the RG equations are given by Uf — U^i = U. As an 
illustrative example, we choose the initial interaction strength {U/A) = 4. 
Our results for this choice of interaction are displayed in Fig. 3. In this 
plot, even though the couplings are initially constant, we demonstrate that 
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they acquire as expected a distinct frequency structure as we integrate out 
degrees of freedom towards the low-energy limit (i.e. A — > 0). Technically 
speaking, this happens because of the interplay of parquet particle-particle 
and particle-hole classes of diagrams in the vertex corrections which are in- 
cluded to infinite order in our field-theoretical RG scheme. In other words, 
despite the irrelevance in a power counting sense of the frequency dependence 
of the coupling functions for the SIAM, they do play an important role in 
the corresponding low-energy dynamics for this case. 

Next we present the results of the RG flow for the quasiparticle weight 
Za of the model as a function of the ratio A/ A. Since at the UV cutoff 
scale (A — Aq) there should be no quantum fluctuations included in the RG 
scheme, the quasiparticle weight must be initially equal to unity (i.e. it must 
assume its noninteracting value) . Naturally, the final value of Za in the low- 
energy limit will depend crucially upon the bare interaction strength U . Our 
results are displayed in Fig. 4 for some choices of initial interaction. In this 
figure, we observe clearly a strong suppression of Za^q with increasing U. 
Despite that, the quasiparticle weight always approach a finite value in the 
low-energy limit. In other words, there is no sign of instability or artificial spin 
symmetry breaking during the RG flow within our scheme. This is consistent 
with the fact that the SIAM conforms to a Fermi liquid description for all 
couplings at low-energies. Our RG approach therefore captures successfully 
this important aspect of the problem. 

In order to verify if the Kondo scale of the model is reproduced by our RG 
scheme, we now plot the inverse quasiparticle weight Z^ 1 in the low-energy 
limit on a logarithmic scale as a function of the dimensionless ratio V jixA. 
This is shown in Fig. 5. For comparison, we also display highly- accurate 
NRG data obtained from Ref. [27] - which yields the correct behavior of this 
quantity for all coupling strengths - and the Bethe ansatz result given by Eq. 
© for the case of (U /it A) > 2 (even though this latter expression is, strictly 




0,0-) ! 1 ! 1 ! 1 1 1 1 1- 

2 4 6 8 10 

A/A 



Fig. 4 (Color online) RG flow of the quasiparticle weight Za as a function of 
A I 'A for several choices of initial interaction strength U within the functional field- 
theoretical renormalization group up to two-loop order. 



12 




U/irA 



Fig. 5 (Color online) Inverse quasiparticle weight Z^ 1 in the limit A — > as 
a function of U/tiA within the functional field-theoretical renormalization group 
(FFTRG) up to two-loop or der . For comparison, we have displayed highly accurate 
NRG data obtained in Ref.|27] and second-order perturbation theory result from 
Eq. (J7J). In addition, we also plot the Bethe ansatz result of Eq. ([5]) for (U/nA) > 2, 
even though this latter expression is only valid in the limit U — s- oo. 

speaking, only valid in the limit U — > oo). Moreover, in order to determine 
the boundaries as a function of the the ratio U/nA which roughly separate 
the weak, the moderate and the strong-coupling regimes of the model, the 
result predicted by conventional second-order perturbation theory (see Eq. 
0) is also added. In this way, we observe in Fig. 5 that since perturbation 
theory results agree well with NRG data for (U /nA) < 1, the value where this 
dimcnsionless ratio is equal to unity could be viewed as marking the boundary 
between the weak and moderate interaction regimes. For (U/nA) > 2.5, 
the exponential (Kondo) behavior of 1/Z emerges, and the strong coupling 
regime of the model finally sets in. From all these remarks, we may conclude 
that our RG results for the inverse of the quasiparticle weight for the model 
are indeed in quantitative agreement with NRG data from weak to slightly 
moderate interactions (i.e., (U/irA) < 1.5), which shows that the present 
implementation of the field-theoretical RG scheme improves somewhat the 
conventional perturbation theory result for this regime (see also the inset of 
Fig. 5). Naturally, it would be also very interesting to test our RG approach 
for the present model in the situation away from particle-hole symmetry, 
where it is well-known that the conventional perturbation theory works less 
well than at the particle-hole symmetric point. Indeed, in the particle-hole 
asymmetric case our RG approach could potentially improve even further 
the conventional perturbation theory result for intermediate interactions. We 
therefore plan to perform this analysis in a subsequent work. 

On the other hand, for (U/nA) > 1.5 in the present case our RG results 
start to deviate quantitatively from the NRG data and the agreement with 
the exact results available for this model becomes only qualitative. This im- 
plies that the exponential suppression of Z contained in the Bethe ansatz 
result is still not reproduced by our RG approach. As can be seen from Fig. 
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5, our results underestimate the strength of fluctuations of the SIAM for 
such large interactions. Therefore, this establishes the limit of validity of our 
two-loop truncation of the RG flow equations for describing this model at 
stronger couplings. Although this approximation proved to be good for not 
too strong interactions, it becomes evidently more and more uncontrolled for 
higher values of U. In order to obtain better results for stronger couplings 
in this model, it would be important to further include higher-order contri- 
butions in the self-energy (e.g., three- loop diagrams or beyond) and also in 
the renormalized coupling functions within the RG truncation scheme. This 
improvement of the results for such large interactions is indeed possible be- 
cause of an important property of the SIAM, i.e. that the system always finds 
itself in a Fermi-liquid state for any coupling strength since it never under- 
goes a phase transition. In view of this fact, a perturbative truncation of the 
functional field-theoretical RG approach about the renormalized impurity 
interaction is expected to converge in this case. 

We can also put our findings into context with some other works available 
in the literature. Karrasch et al. recently implemented a Wilsonian functional 
RG study [35] of some quantities of the SIAM such as the effective mass, 
the static spin susceptibility and the spectral function. They concluded that 
while their RG method also reproduced precise NRG results for these phys- 
ical quantities at weak couplings, they failed to capture the correct Kondo 
scale generation for the model within their truncation. In another recent 
work [28], a functional RG scheme within the Keldysh formalism was im- 
plemented, and the authors obtained quantitative agreement with accurate 
NRG data for the SIAM in the equilibrium situation only for weak inter- 
actions. All these findings are in good agreement with our results for this 
regime. In addition to this, since the functional field-theoretical RG method 
described here turns out to be easier to implement at higher loops [111 than 
the Wilsonian approach, higher-order calculations within the present scheme 
are indeed feasible and could improve the results for this model at stronger 
couplings. In view of this, we believe the present RG scheme can be viewed as 
a possible alternative to other functional RG methods to describe equilibrium 
properties of the SIAM. 

5 Conclusion 

We have implemented a functional field-theoretical RG method up to two- 
loop order to the symmetric SIAM at zero temperature within the wide-band 
limit. One of the objectives of this study was to extend the conventional field- 
theoretical RG scheme to include the dynamics of all frequency-dependent 
renormalized parameters of the model. We have shown that, in order to 
accomplish this, it is convenient to define a soft ultraviolet regulator in the 
space of Matsubara frequencies for the renormalized Green's function of the 
model. In this way, we have derived analytically and then solved numerically 
the corresponding RG equations both for the coupling functions of the model 
and the quasiparticle weight in the low-energy limit. These equations fully 
included the interplay of particle-particle and particle-hole parquet diagrams 
and also the effect of two-loop self-energy feedback into them. Naturally, 
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the RG framework was ideally suited for achieving this goal in an unbiased 
manner. 

One important point we want to stress here is that, because of the two- 
loop self-energy contribution into the RG flow equations, our approach did 
not suffer from any instability (or artificial spin symmetry breaking) for any 
choice of interaction strength U. This is consistent with the fact that the 
SIAM conforms to a Fermi liquid description for all couplings at low-energies. 
Our RG approach therefore captured successfully this aspect of the problem. 
In order to assess the quality of our two-loop truncation of the RG flow 
equations, we have benchmarked our results against highly accurate NRG 
data obtained from Ref. . As a result of this comparison, we have con- 
firmed that our results indeed reproduce quantitatively the NRG data for 
the model from weak to slightly moderate interactions. Moreover, since the 
functional generalization of the field-theoretical RG method described here 
turns out to be technically easier to calculate at higher loops as compared 
with the more conventional Wilsonian RG approach, higher-order computa- 
tions within the present scheme are relatively straightforward to implement. 
This latter calculation could potentially improve our results for this model 
at stronger couplings. Therefore, we believe the development of the present 
functional field-theoretical RG approach could offer a promising alternative 
to other functional RG methods which aim to describe electronic correlations 
within the Anderson impurity model. 
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A Appendix 



In this appendix, we write down the explicit form of all the functions defined in 
the text. They are given by 



A-f (U>1,U>2,W3) 



r 



(w + Z/iZ\sgn(cj))(oji + L02 — oj + ^A^4sgn(a;i + L02 — w)) 



(ui,U2,U)), 



(17) 



_B^(cJl,CJ2,OJ3) 



[ 



dw- 



-T 2 (u)i,u) 3 ,u) 



(w + ZAAsgn(u))(uji — CJ3 + u + ZAAsgn(uji — us + w))* 

U^(U1,U,U33)U^(U1 — OJ3 + U, LO2, Ul) + Ufi(ui,U),(d3)Ufi{ui —(jJs+LU, LJ 2 , w)J , 

(18) 
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/oo ^ 
.oo (w + Z^^sgn(w))(w + W2 - ^3 + Z/i/Asgn(u; + cu 2 - w 3 )) 

[i/-fj.(wi, W + CJ 2 - W3, W)?7ti(w, W2, W3) + U-f (Wl , W + W 2 - W3, w)i/-f 1 (w, W2, W3) 

— Uf(ui,U + LU2 — U)' A ,U))U^(u), U)3,L0 + L02 — LO3 ) — (cOl , CO + C02 — OO3 , UJa)U^ (<jJ , 0J2, OJ3) 

+Z7f4,(wi,W + ^2 - W3,w){/-f4,(W2,W,W3) + (Wl , W + W 2 - ^3 , W4) ( W > ^ , W3)] , (19) 

/oo ^ 
dti; 7 — : 7 /a — r^v — ; 1 ? a — i — ; v\- 77 i 1 ( Wl ' W2 ' w ) 
_oo (w + Z^/lsgn(w))(wi +L02 - co + Z A Asgn(ui + lo 2 - to)) 

[{/•A(W1, W2, W){/A( W 1 + ^2 - W, W, W3) + ^A( Wl > W2 > W )^A( W ' Wl + ^2 - W, W4)] , 



(20) 



/oo 
(w + Z A Asgn(u))(coi - U3 + w + Z^Z\sgn(o;i - o; 3 + w)) 
[t^( Wl ' w ' ^3)^(^1 — W3 +w,cj2,^) + ^(cji,cj,u;3)^t4,(^i - 0J3 + w,u;2, w)j , 



(21) 



dtJ ? — 1 7 a — r^v — ; 1 7 /v — ? — ; vT-^s 1 ^, w 3 ,w) 

(w + ZAAsga[u))(u + U2 - w 3 + Z/iZlsgn(cj +w 2 - W3)) 

[{/A(W1, W + ^2 - W3, W4)C^A( W ' W2 > W 3) + C/A(^l, W + ^2 - W3, w){7f4.( W > W2 > W + W2 - Ws ) 
-[/A(W1,W + W2 - L03,L0)U^(U!,C02,0J3)^ , (22) 



where 



7?(wi,W 2 ,w) = f [M + |W1 +W2 -Wn/fdwl/yl^dwi +W2 - W|/A), 

^(wi.ws.w) = f [M + |wi -uj3+oo\]K(\oj\/A)K{\co 1 -lo 3 + co\/A), 

J^(W2,W 3 ,W) = fj^) [M + |W+W2 -W 3 |]«'(M/A)X'(|W+W2 -W 3 |M). 



(23) 
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